Metal-insulator transition in two-dimensional disordered systems with power-law transfer terms 
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We investigate a disordered two-dimensional lattice model for noninteracting electrons with long-range 
power-law transfer terms and apply the method of level statistics for the calculation of the critical properties. 
The eigenvalues used are obtained numerically by direct diagonalization. We find a metal-insulator transition 
for a system with orthogonal symmetry. The exponent governing the divergence of the correlation length at the 
transition is extracted from a finite size scaling analysis and found to be = 2.6 ± 0. 15. The critical eigenstates 
are also analyzed and the distribution of the generalized multifractal dimensions is extrapolated. 
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Starting with Anderson's worklll more than four decades 
ago, the disorder driven metal-insulator transition carries on 
to be a subject of active research. The underlying basic con- 
cept has proven to be important in many physical situations 
because it provides an essential requirement for the micro- 
scopical understanding of certain quantum phase transitions 
in condensed matter physics. According to the scaling the- 
ory of localization^ noninteracting electrons in infinite one- 
dimensional (ID) and two-dimensional (2D) disordered sys- 
tems are always localized at zero temperature. It is known that 
there are exceptions to this rule. In particular, the quantum 
Hall system (absence of time reversal symmetry) with its di=. 
verging localization length in the center of the Landau band* 
and the 2D disordered system with symplectic symmetrycE 
are notable examples. However, even in the presence of time 
reversal symmetry, indications of a metal-insulator transition 
have been reported inpcxperiments on two-dimensional elec- 
tron and hole systems.Q These findings have further intensified 
the efforts to look for possible mechanisms that allow a metal- 
insulator transition in systems with spatial dimension d < 3. 
Clearly, electron-electron interactions that were not taken into 
account in Ref. 2 are important at low temperatures and pos- 
sibly may account for the observed behaviorQ 

Another class of systems exhibiting a metal-insulator tran- 
sition for noninteracting particles in dimensions d < 3, can be 
represented by models that in contrast to the Anderson model 
include long-range transfer terms. |-This possibility was al- 
ready indicated in Anderson's paperjllbut later addressed only 
a few times. For instance, the decay of the localized eigen- 
states has been investigated numerically for a ID tight-binding 
system with long-range transfer termsE Also, the localiza- 
tion and dephasing of dipolar excitons in topologically disor- 
dered systems and the influence of weak longijguige hopping 
on the 3D Anderson model have been studiedBEiM as well as 
the problem of Bloch states for a particle moving fast through 
a lattice of Coulomb scatterers with power-law singularity.llil 
Instead of triggering the transition by varying the disorder 
strength, the crossover from localized to extended states has to 
be tuned by changing the pa^pj^npnt (3 of the power law transfer 
terms Vi 
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Meanwhile, it has become clear that the behavior at the crit- 
ical point and its multifractal eigenfunction statistics vcmi be 
described by a one-parameter random matrix ensemble.Ej An 
intensively studied example is the power-law random banded 



matrix model (PRBM),i3l^ which exhibits a transition from 
localized to extended states and contains multifractal eigen- 
states at the critical point. This was demonstrated by map- 
ping the PRBM onto a nonlinear sigma model, which could 
be solved for limiting cases. In one thereof, the PRBM trajo^ 
forms to the ID version of the model studied by Levitov.Erty 
The corresponding critical properties have been investigated 
numerically recently .0 

Despite this considerable success, up to now it was not pos- 
sible to solve the experimentally more important 2D mod- 
els with long-range transfer terms analytically. Therefore, no 
quantitative estimates have been obtained for the critical quan- 
tities like the exponent v of the correlation length or the fractal 
correlation dimension D{2). The latter essentiaJly influences 
the dynamics at the metal-insulator transitionllJ In this paper 
we present the first results for the multifractal properties and 
the critical exponent of a 2D system with orthogonal symme- 
try, which was found to be = 2.6 ± 0.15. 

The model investigated here describes noninteracting par- 
ticles on a 2D square lattice with distance dependent random 
transfer terms Vij = Vsij {b/\ri — rj\)^ and random on-site 
potentials e^. The corresponding Hamilton operator is given 
by 



ri){ri 
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The {si} and {e^ } are two sets of uncorrected random num- 
bers uniformly distributed between — W/2 and W/2, and be- 
tween -5/2 and S/2, respectively. We take b ^ 1, S/V 1, 
and fix the diagonal disorder by W/V — 6. The unit of en- 
ergy is set by y = 1. The special choice of the diagonal 
disorder strength W places the system in a certain volume of 
the parameter space where the transition is easier to access. 
The transfer terms, which allow transitions from a given lat- 
tice site to all others, depend on the corresponding distance 
via the power-law \ri — rj\^^. By changing j3 one tunes the 
metal-insulator transition. L is the size of the square system 
measured in lattice spacings a, and Dirichlet boundary condi- 
tions are applied in both directions. 

To investigate the localization behavior we apply the level 
■etbpd which has proven to be a very powerful tool 
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in the past.t^lfJ The required eigenvalues for the level statistics 
and the eigenvectors for the multifractal analysis were calcu- 
lated by direct diagonalization. For the above chosen ratio of 
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the diagonal and nondiagonal disorder, the density of states is 
symmetrical about the band center E /V = and only weakly 
energy dependent within the range [—2V,2V]. The number 
of realizations calculated was such that for each set of param- 
eters the total number of eigenvalues on the average add up to 
6 • 10^ for system size L/a < 40 and 3 • 10^ for larger sizes. 
As in previous work, we choose the quantity Io{f3, L) as our 
scaling variable because this, choice does not depend on any 
arbitrary cutoff parameter£0 Here, /o(/3, L) = (s^)/2 is half 
the second moment, (s^) = s^Pp j^ds, of the probability 
density P[}.l{s) for finding an energy difference of two con- 
secutive eigenvalues, s — \En — -En+il/A, where A denotes 
the mean level spacing. We performed a proper unfolding pro- 
cedure and checked that the results were independent of the 
width of the energy interval around E/V = Q from which the 
eigenvalues were taken. 

The eigenvalue statistics has been calculated within the en- 
ergy interval [— l.SV, 1.5^] for square systems of linear size 
L/a = 10, 15, 20, 25, 30, 40, 45, 50, 55, 60, 70, 80, 100, and 
120. The exponent (3 of the power-law decay of the trans- 
fer terms was varied in the range 1.7 < (3 < 2.3. A total of 
284 Jq L) values have been accumulated. Increasing the 
size L, the magnitude of Iq decreases for /? < 2, but gets 
larger for (3 > 2 with an almost scale independent value at 
the point of intersection near (3 = 2. The latter behavior is a 
signature of a critical level statistics connected with a quan- 
tum critical poinnlj. In the limit L oo, the scaling variable 
Iq can take on three values. For j3 < 2 only extended states 
are expected so that /q = 0.643 for orthogonal symmetry that 
is the universal random matrix result for the diffusive regime. 
For /3 > 2 all states will be localized in the thermodynamic 
limit so that Iq = 1.0 which corresponds to the Poisson proba- 
bility density distribution, P{s) = exp(— s). This means that 
the probability of having two neighboring eigenvalues close 
together is maximal, which is due to the negligible overlap of 
the associated localized eigenvectors. The third possible value 
Iq can assume in the limit L ^ oo is the critical value Iq at 
/3c = 2, which is nonuniversal and known to depend on model 
parameters as well as,on the boundary conditions and on the 
shape of the system.E3 

In Fig. |l], which for clarity shows only a subset of the raw 
data, a crossing of the /o(/9; L) values belonging to different 
system sizes L can be recognized near (3 = 2. However, the 
weak monotonous shift of the crossing point with increasing 
system size L indicates an irrelevant scaling term. From the 
Io{P = 2, L) data we get an estimate for the irrelevant scal- 
ing exponent y (see inset of Fig. [l|) which then was finally 
determined together with the correlation exponent i/ using a 
renormalization group ansatz where the_dimensionless quan- 
tity Ioi(3, L) is expressed as a function,El 



IoW,L) = MA.L'^n + A2LyhiA,L'^n, (2) 



containing a relevant and the leading irrelevant scaling vari- 
able Ai and A2, respectively. The corresponding scaling ex- 
ponents are ly and y. The functions /i, /2 as well as the scaling 
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FIG. 1: /o(/3, L) of the level statistics for a 2D disordered system 
with power-law transfer terms. The system size L/a is 10 (•), 
15 (x), 20 (*), 25 (□), 30 (+), 40 (o), 60 (■) for power-law ex- 
ponents in the range 1.7 < /3 < 2.3, and L/a = 80 (A), 100 (O) 
with 1.9 < (3 < 2.1, and 120 (V) with 1.925 < (3 < 2.075. The 
inset shows the size dependence of the /o(/3 = 2, L) data which can 
be fitted by the relation Io{(3 = 2, L) = IS + cL^ with /J ~ 0.78 
and y ~ —1.2. 



variables are expanded in Taylor series with 7 = (/? — (3c) / (3c, 



Io{(3,L) 2, J2iCp + LyDp){jL'/n^ (3) 

Here, Cp and Dp are the coefficients that control the relevant 
and irrelevant scaling fields, where Co = Iq determines the 
critical value of Io{(3, L) at (3c. Assuming that the deviations 
between the data and the model are random and caused by un- 
correlated statistical noise due to the inherent disorder and the 
finite number of realizations, the best fit was attributed to the 
set of parameters with minimaLi^ statistics. Since the critical 
power-law exponent is known,EHlB (3c = d, we can fix (3c = 2 
and obtain for the orthogonal 2D case /g = 0.784 ± 0.005, 
and a critical correlation exponent ly = 2.6 ±0.15. The 
leading irrelevant scaling term is governed by an exponent 
y = — 1.25 ± 0.25. The resulting scaling curve is shown 
in Fig. ^ where the original data corrected for the irrelevant 
scaling terms, 7o(/3, L) = h{(3, L) - ^^=0 ^5^(7^'/")^ 
are plotted versus \(3 — 2\L^^^ . We achieve almost the same 
results if only 178 data points from the restricted interval 
1.8 < /3 < 2.2 are taken into account. In this case the 
quadratic terms in Eq. (S) can be neglected. To test whether 
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FIG. 2: Scaling of the variance Iq{[3,L) of the level statistics for 
a 2D disordered system with power-law transfer terms. The system 
size L/a varies between 10 and 120 and the power-law exponent is 
taken from the range 1.7 < P < 2.3. The finite size scaling analysis 
yields Jq = 0.784 ± 0.005, a critical exponent u = 2.&± 0.15, and 
an irrelevant scaling exponent y = —1.25 ± 0.25. 



our ansatz for the correlation length, £, ^ \I3 — /5c |^^, is justi- 
fied at all, we directly calculated ^{0) from our adjusted data 
using the method described in Ref. 24. The result, which is 
shown in Fig. ||, is consistent with our assumption and the best 
fit to the calculated dala yields /i = 2.4±0.1. It is well known 
that with this methods one gets an estimate only for the lower 
bound of the true critical exponent v. It is important to men- 
tion that our result for the correlation length in 2D differs 
completely from the asyip»ietrical function proposed for the 
corresponding ID model.ll3lJ 

Our final result for the critical exponent, v = 2.6 ± 0.15, 
is larger than the value ^accepted for the quantum Hall sys- 
tem u = 2.35 ± 0.03,0 but lies within the range of those 
reported for systems with symplectic symmetry (see, e.g., 
Ref. 25). Also, we would like to remark that our result for 
the correlation exponent is in good agreement with the value 
zv — 2.6 ± 0.8 obtained recently from temperature scaling 
in an experiment on a 2D system where the disorder was gen- 
erated by a layer of self-assembled quantum dots.Ej Thereby, 
one has to assume z = 1 for the dynamical exponent as ex- 
pected for interacting electrons. However, there are in the lit- 
erature also reports of smalletcJ as well as larger valuesEa for 
the coiTelation exponent v in 2D electron and hole gases. Pra=. 
vided that the accordance between the experimental valueEl 
and our result for the model presented above is not fortuitous, 
the coincidence indicates that both metal-insulator transitions 
belong to the same universality class. 

Now, having determined the critical exponent of the cor- 
relation length for our model with power-law transfer terms. 




FIG. 3: The dependence of the correlation length £,{13) on the power- 
law exponent /?. The values are, up to a prefactor, calculated directly 
from the Io{P, L) data which are the raw data /o(/3, L) adjusted for 
irrelevant scaling terms. The dashed line shows the best fit f ~ 



we turn to the analysis of the critical eigenstates (/3 = 
2) for a square system of linear size L/a = 150 and 
eigenenergy near zero (band center). A multifractal analy- 
sis has been performed along the standard procedure where 
the scaling of a "box-probability" is calculated, P{q, A) — 

Ef ^'HE.eo.(/) IV'£(r)n' ~ A^('), from which the gen- 
eralized fractal dimensions D{q) — T{q)/{q — 1) or, by a 
Legendre transform, the so called /(a((7))-distribution can 
be derivedo Here, is the ith box of size I = XL 

from which the qth moment of the modulus of the normal- 
ized eigenstate i/jEir) is taken. In Fig. ^we show the f{a{q)) 
distribution in comparison with the parabolic approximation, 
f{a{q)) =d~ {a{q) - a(0)) V(4(ci(0) - d)), which for fi- 
nite systems is valid only for small |g|.E] As usual in studies on 
finite systems, the deviations are larger for such a{q) that cor- 
respond to negative q values because the negative exponents 
blow up those spots where the eigenstate almost vanishes. The 
correlation dimension is extracted to be D{2) = 0.9 ± 0.05 
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FIG. 4: The /(Q(g))-distribution of a critical eigenstate with energy 
at the band center for a square system of linear size L/a — 150. The 
parabolic approximation f{a) = d — {a — a(0))^/ (4(a(0) — d)) is 
fitted by a single parameter a(0) = 2.92 ± 0.05. 
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so that the exponent rj = d — D{2) ~ 1.1 describing anoma- 
lous diffusion near the critical point is very large compared 
to the other known 2D situations, e.g., rj ~ 0.38 for the 
quantujn Hall systerro and -q ~ 0.35 for the symplectic 
case.tll Therefore, one has to expect a strong influence on the 
dynamical properties at the metal-insulator transition origi- 
nating from th^patial amplitude fluctuations of the critical 
eigenfunctions.El 

In conclusion, we investigated a disordered two-dimen- 
sional lattice model with long-range transfer terms which are 
governed by a power-law decay ^ jr^ — rj\~^ . Using the 



method of level statistics, the existence of a metal-insulator 
transition has been demonstrated in a 2D system with orthogo- 
nal symmetry. The transition is tuned by the power-law expo- 
nent /3 of the transfer terms. From a finite size scaling analysis 
the critical exponent determining the divergence of the corre- 
lation length near the transition has been obtained. We found 
V ~ 2.6 ±0.15 which is very close (assuming a dynamical 
exponent z = 1) to an experimental result, zv — 2.6 itjO.8, 
which was estimated recently from temperature scaling.ES 

We thank Ferdinand Evers for very helpful discussions. 
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